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ABSTRACT 

We report the first detection of the ground-state rotational transition of the methylidyne cation CH + towards the massive star-forming region 
DR 21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad 
absorption features associated with the DR 21 molecular ridge and foreground gas. These observations allow us to determine a 12 CH + y=l-0 line 
frequency of v=835137+3 MHz, in good agreement with a recent experimental determination. We estimate the CH + column density to be a few 

00 ' 10 13 cm" 2 in the gas seen in emission, and > 10 14 cm" 2 in the components responsible for the absorption, which is indicative of a high line of sight 
average abundance [CH + ]/[H]> 1.2 x 10" 8 . We show that the CH + column densities agree well with the predictions of state-of-the-art C-shock 

' . ' , models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines. 
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1. Introduction 

The methylidyne ion CH + was among the first molecules to be 
detect ed in the interstellar medium (ISM) Douglas & Herzberg 
(Il94lh . This reactive ion is prevalent in the diffuse ISM with 
column densities several orders of magnitude above the pre- 
dicti ons of U V-driven equilibrium models (see references in 
Godard et al. 2009). Apart from the strong emission lines of 
CH + detect e d in both the envel ope of the Red Rect angle 
iHobbs et all J200l and NGC 7027 lCernicharoet alJ Jl997b . all 
other observations are absorption lines detected at visible wave- 
lengths in the spectra of nearby stars. CH + is a light molecule 
and its rotational lines lie at submillimeter and far infrared wave- 
lengths. The exact frequencies of the rotational transitions of 
CH + have remained elusive for a long time because of the ex- 
treme reac tivity of CH + and the difficulty in isolating it in the 
laboratory iPearson & D rouin (2006). Recent laboratory mea- 
surements have led to v = 835.137498(20) GHz for the ground- 
state transition (Amano, 2010). Ground-based astronomical de- 
tection of 12 CH + (l-0) is prevented by its proximity to a strong 
atmospheric line of water vapor. The ground-state frequency of 
the isotopologue 13 CH + , redshifted by ~ 5 GHz, has superior sky 
transmission. It has been detected at the Caltech Submillimeter 
Observatory towards several massive star - forming regions of th e 
inner Galaxy (Falgarone et al. in prep., iFalgarone et al.l l2005). 
The 12 CH + abundances averaged along these long lines-of-sight 
(los), confirm the high abundances of this species inferred from 
visible observations in the local ISM, [CH + ]/[H]~ 8xl0" 9 on av- 
erage. In this paper, we report the detection of the 12 CH + 7=1-0 
transition towards the massive star- forming region DR21, pre- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



sented in Sect. 2. The HIFI observations are described in Sect. 3. 
The results, given in Sect. 4, are compared to models in Sect. 5. 

2. The DR 21 region 

The massive star- forming region DR 21 is located in the Cygnus 
X complex at an average distance of 1 .7 kpc, that of the 
Cyg OB2 stellar association (see reviews by iJakob et alj|2007l: 
Schneide r et al.1 t2010l) . The DR 21 molecular core is one of 
the most mas sive cores in the Galaxy (2.5xl0 4 M at the 
1 pc-scale) (Kirby 2009). It is located in front of the main 
DR 2 1 HII region co mprising five compact HII regions and more 
diffuse ionised gas (Ro elfsemaet al.1 1 19891: 1 Wilson et all 1 19951: 
Cvganowski et al. 2003). 

The star- forming region DR 21 itself is known to host one of 
the most power ful molecular outfl ows shining in vibrationally 
excited H2 (e.g., Garden & Carlstrom 1992). The outflow source 
is heavily extinguished by more than 10 magnitudes, and its lo- 
cation is not accurately known dWilson et al.ll 19951) . High veloc- 
ity HI (up to 90 kms" 1 ) associa ted with the molecu lar outflow 
has been detected by the VLA dRusseil et al.lll992l) . indicating 
that the molecular outflow may be driven by this atomic je t. The 
outflow shows up as broad wings in CO and HCO + lines iKirbvl 
(2009). At the position of the HIFI beam, the v = 1-0 S(l") 
H? em i ssion line has a simi l ar bro ad width iNadeau & Belandl 
( 1988); C ruz-Gonzalez et al.l d2007l) to that of the SiO emission 
( Mott e et al.ll2007l) . 

Some line profiles imply that there is foreground gas, which 
is probably associated with th e Cygnus X complex IJakob et all 
(120071) : [Schneider et alJ d2010h . A foreground component asso- 
ciated with W75N shows up as weak emission lines at v^fi ~ 
8 km s" 1 in low excitation transitions of CO and HCO + or 
atomic carbon, and in t he atomic oxygen fine- structure line, with 
N(Q) > 5 x 10 18 cm" 2 iPoglitsch et all (Il996l) . The total hydro- 
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Fig.l. [CII], CH + (l-0) and p-H 2 Herschel/HIFI spectra (top 
panels) and HCO + (1-0 ) IRAM-30m spectrum (bottom panel) 
(ISchneider et"ai] l2010h . The HIFI spectra are taken in DSB, 
hence saturated lines reach to half the continuum level (dotted 
lines). The weak line in the C + panel is that of an intermediate 
OFF position. The CH + (J=l-0) spectrum is shown assuming 
a rest frequency of 835137 MHz. The Gaussian velocity com- 
ponents involved in the empirical model of the CH + profile are 
shown (dotted blue curves) as well as the resulting models for 
all lines (red curves) (see Sect. 4). 



gen column density of this foreground material is estimated to 
be N{H) ~ 1.4 x 10 2 2 cm" 2 for an eleme ntal abundance ratio 
[0]/[H] =3.45 x 10" 4 lOliveira et all (120051) . in coherence with 
N(U) - (1.3 + 0.4) X 10 22 c m' 2 inferred from K extinction to 
DR21 ((Marshall et a l. 2006). This component is also detected 
in atomic hydrogen, which indicates t hat there is sa turated ab- 
sorption of between -5 and 18 kms~' Rob erts et al.l d 19971) . Its 
estimated HI column density, if T s = 2QK, is N(HT) > 1.5 x 10 21 



-2 



cm 



3. HIFI observations 

All FIR spectra presented here were obtained in the course of 
performance verification observations with the HIFI instrument 
(Pilb ratt et alJ 120101 : Ide Graauw et al. 2010). Since their main 
goal was to demonstrate the functionality and performance of 
the different observing modes, the spectra were taken with a 
large variety of observing strategies. Most observations were 
only single-point observations towards the central position of 
the DR21 HII region at RA=20h39m01.1s, Dec=42°19'43.0" 
(J2000). At the frequencies around 1 THz discussed here, the 
Herschel beam covers about 20" HPBW (or 0.16 pc at 1.7 kpc). 
Here, we only present data from the wideband spectrometer 
(WBS) that provides a resolution of 1.1 MHz, corresponding to 
0.2 km/s (at 1900 GHz) and 0.4 km/s (at 835 GHz). 

The CH + spectrum at 835 GHz was obtained in a single- 
point load-chop observing mode without OFF reference. It uses 
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Fig. 2. Superimposition of the CH + (red), C + (black), and p-H 2 
(blue) spectra (from which the continuum has been subtracted) 
to illustrate the good coincidence of the blue wing of the line 
profile when the CH + (l-0) frequency is v = 835.137 GHz. We 
note also the coincidence of the absorption features in CH + (l-0) 
and C + . 

the internal cold load as reference, having the advantage that we 
can exclude any self-chopping effects in the line profiles. The 
mode has, however, the disadvantage that the different optical 
paths towards the sky and the internal calibration load sometimes 
lead to standing wave differences that are detected as baseline 
ripples. We subtracted these manually using the HifiFitFringe 
pipeline tool. The total integration time on the source was 48 s 
leading to a noise level of 0. 1 K for the combination of the spec- 
tra from both polarisations at the WBS resolution. 

The H 2 1 1 1 -Ooo (1113 GHz) line was observed in the frame 
of a spectral scan in a load-chop mode and an on-source integra- 
tion time of 30 s. The baseline was calibrated on a distant OFF 
position (20h37ml0s, 42°37'00") because the first selected OFF 
position, closer to the ON, retained some signal (see Fig. 1). The 
[CII] line was detected in one of the first double-beam switch 
raster mapping observations in which neither the logic of the 
mapping mode worked as expected nor the pointing of the in- 
strument was yet known. The instrument pointing was measured 
afterwards and the map coordinates were corrected by hand to 
include the correct offsets. A comparison of the integrated [CII] 
intensities in the shifted map with the MSX 8 jt/m band shows a 
very good agreement, confirming the accuracy of the correction. 
The observations were taken on a 7" raster that included a 14 s 
on-source integration time at each point. To obtain a spatial res- 
olution comparable to that of the other two lines, the map was 
convolved with a Gaussian beam of 20" at the DR 21 central 
position. 

4. Results of the line profile analysis 

The HIFI spectra are displayed in Fig. [TJ with, for comparison , 
an IRAM-30m spectrum of HCO + (l-0) lSchneider et alfeOlOh . 
The line profiles exhibit both complex emission and absorption 
because of the multiplicity of velocity components and gas phys- 
ical conditions within the beam. The line integrated areas in 
emission and absorption are remarkably comparable. This would 
lead to a weak or non-detection with low spectral resolution in- 
struments. The dashed lines in the CH + and P-H2O spectra, at 
about half the continuum level, reveal the broad velocity ranges 
across which the absorption lines are saturated. 

The CH + (l-0), C + and p-H 2 spectra, plotted after their con- 
tinuum has been subtracted, are superimposed in Fig. [2] for the 
frequency of the CH + /=l-0 line set to 835.137 GHz. We pro- 
pose in Sect. 5 that the broad line emission originates from a 
shock associated with the outflow, and we show that the three 
species, CH + , C + , and H2O are connected with each other in the 
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Table 1. Empirical models of the broad emission profiles 



Table 2. CH + Column densities 



Line 



kms 



Av 
kms -1 



T 

K 



CH+ 


-3.0 


20.0 


2.2 


p-H 2 


-1.0 


20.0 


2.16 


C + 


-1.8 


20.0 


22 







80.0 


1.0 


HCO+ 


-3.6 


17.0 


2.8 



shock chemistry. The good agreement (within + 1 kms 1 ) be- 
tween the blue wing of the emission line profiles of these three 
species allows us to adopt v = 835137 + 3 MHz as the rest fre- 
quency of the 12 CH + (l-0) line, which is remarkably consistent 
with the experimental value measured by Amano (2010). 

The CH + 7=1-0 line profile can be decomposed into a broad 
emission line and a series of absorption features. We used infor- 
mation from the line profiles of CO and HCO + 7=1- (see Fig. 
[TJ, the spectr a of atomic car bon Jako b et all (120071). and ot her 
HIFI spectra (lOssenkopf et al.ll2010t Ivan der Tak et alj|2010l) to 
constrain the velocity and linewidth of these components. Since 
the I3 CH + line has not been observed, the actual CH + emission 
profile is unknown. We had to rely on the symmetric shapes of 
the CO emission lines to model the broad emission wings with 
Gaussians adjusted on the blue wings of the lines (Table 1). Over 
the velocity range -5 to 17 km s , not only the continuum level 
but also the core line emission is absorbed by intervening gas. 
The broad emission lines (Av = 20 kms -1 , see Table 1) cor- 
respond to the emission associated with the outflow shock (see 
Sect. 5). A still broader component in the [CII1 line has a width 
similar to that of the H66ff recombination lines lRoelfsema et al.l 
(1989). Self-absorption features in the range ~ -3 to 2 kms -1 , 
caused by the close environment of the DR21 core, are visible 
in CO, HCO + (4-3), and other dense gas tracers (e.g., NH3 inver- 
sion lines). The absorption dip at this velocity is more prominent 
in the H20(lij - Oo,o) profile as expected from the high critical 
density of this transition. We disregarded this velocity range in 
our analysis because of the unknown line intensity of the dense 
core emission. In contrast, the absorption at v > 7 kms -1 cor- 
responds to the gas that is most likely to be associated with the 
W75N cloud in the Cygnus X complex. Its velocity coverage is 
very similar in the CH + and C + spectra. We fitted three compo- 
nents in this velocity range, guided by the absorption features of 
the HCO + 7=1-0 and H z O profiles at v LSR = 7.5, 11, and 13.7 
km s -1 shown in Fig. 1, the resulting profiles being in red. 

The CH + column density of the emission line was deter- 
mined by assuming that it is optically thin. This is consistent 
with the 13 CH + (l-0) line not being detected in emission towards 
galactic massive star-forming regions (Falgarone et al. in prep., 
Falgarone et al. 2005). The total CH + column density was there- 
fore assumed to scale with the line integrated area (in Kkm s -1 ) 
with a coefficient that depends on the excitation temperature. For 
40 K< T ex < 100 K, the range of T ex that minimizes N(CH + ), 
a lower limit N em (CH + ) = 7 x 10 n cm -2 JT(v)dv was ob- 
tained, using the spontaneous decay rate A 10 = 5.9 x 10~ 3 s -1 
inferre d from the dipole moment, \i = 1 .62 D IKobavashi et al.l 
(Il993l) . In absorption, we inferred the column density from the 
integral of the optical depth (in kms -1 ) to be Af a fc(CH + ) = 
3.3 x 10 12 cm -2 J T(v)dv. It is almost independent of T ex as long 
as T ex « hv/k = 39K, an approximation that is valid for gas 
of density much lower than the critical density of the transition 
(> 10 6 cm -3 ) and not closely associated with intense FIR radi- 
ation. The resulting CH + column densities are given in Table 2. 



Component 



K 



N(CH + ) 



Broad line emission 
Foreground absorption 



40K < r„ < 
3.0 



100K 



2.2 x 10" 
> 1.7 x 10 14 



Those of the saturated absorption components are lower limits. 
Although the Gaussian decomposition is by no means unique, 
it provides an estimate of both the detected column densities 
and the CH + abundance relative to the total hydrogen [CH + ]/[H] 
> 1 .2 x 10 -8 in the diffuse gas. We note that the [CII] line opacity 
of the foreground gas, averaged over 10 kms -1 , is Tfcrrj > 1.8. 
This opacity limit corresponds to a C + column density larger 
than 2 x 10 1S cm -2 for gas of density lower than the critical den- 
sity of the transition w cf (C + ) , i.e. ~ a few 10 3 cm -3 , whatever the 
line excitation temperature (Crawford et al. 1985). The required 
C + column density is much larger if the transition is thermalized, 
i.e. n » n cr {C + ). From the estimated gas column of 1.4 x 10 22 
cm -2 , the expected C + column is ~ 2 x 10 18 cm -2 , assuming a 
gas-phase carbon abundance of 1 .4 x 10 -4 . This value is in close 
agreement with that derived for the low density limit, strength- 
ening the association of the foreground absorption with gas of 
densities lower than ~ a few 10 3 cm -3 . 



5. Comparison with model predictions 

The only formation route of the molecular ion CH + is under- 
stood to be initiated by the reaction C + + H2 — * CH + + H, 
which is highly endoenergic (AE/k =4640 K). The formation of 
CH + in the cold interstellar medium therefore requires suprather- 
mal energy. Several scenarios hav e been investigated : C-shocks 
dPineau des Forets et all 1 986 (@), highly illuminated and dense 
photon-dominated r egions (PDR) where C + reacts with vibra- 
tional ly excited H2 dSternberg & D algarno 1995; Ag undez et al.l 
2010), tu rbulent interfaces b etween the warm and cold neutral 
medium dLesaffre et alJ l2007). an d regions of intermittent turbu- 
lent dissipation (TDR models. iGodard et al.ll2009t) . In the first 
and last models, ion-neutral friction plays a major role. 

For the DR 21 core, it is unclear which processes are at 
work in the CH + formation. Broad (Av = 20 kms -1 ) emis- 
sio n lines of S iO Motte et al. (2007) and vibrationally excited 
H? ICruz-G onzalez et al. (2007), similar to the broad emission 
in the HIFI profiles, are observed and suggest the influence of 
a shock. We now present a preliminary attempt to account for 
the observed column densities of CH + in this highly irradiated 
environment. 



5.1. CH + emission line: C-shock models 

We used a steady-state mod el for C-shocks with time-dependen t 
chemistry and ionization dFlower & Pineau des Foretsl I2003). 
The chemical network was supplemented with the relevant 
photo-reactions (Pineau des Forets, priv. comm.). The H2 and 
CO self-shielding were assumed to be those of a PDR model at 
«h = 10 4 cm -3 . We adopted a pre-shock density n# = 10 4 cm -3 , 
a shock velocity v s =20 kms -1 , a ma gnetic field B = 2 00jLfG 
(from the measure of magnetic field in Ro berts et al.l fl997). and 
a standard cosmic-ray ionisation rate of f = 5 x 10 -17 s -1 . We 
stopped the computation of the shocks at a neutral flow time of 
10 4 yr, which corresponds to the estimated dynamical age of the 
outflow (lobes of 0.5 pc for a jet velocity up to 100 kms -1 ). 
The CH + column-density was found to have already reached a 
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Fig. 3. CH + (blue) and H2O (red) column densities integrated 
through the C-shock model (see text) for a broad range of UV- 
field irradiations^-. The shielding of the shocked material is A,, = 
1 . We note the opposite behaviours of CH + and H2O with x- 



plateau at this stage in the shock, so the CH + results are insensi- 
tive to that choice of age. 

Although the incident radiation field is estimated to be^f = 
10 5 in ISRF units lOssenkopf et al.ld2010l) for the material closest 
to the HII region, we keep this value as a free parameter. Figure 
[3]displays the total column densities of CH + and H2O as a func- 
tion of the adopted external radiation field. Its effect turns out to 
be quite interesting: a strong radiation field (x > 10 3 ) is needed 
to account for the large observed CH + column density. It is the 
enhanced photo-dissociation of CHj and CHJ , both products of 
the rapid hydrogenation of CH + in the shock, that increases the 
abundance of CH + as x increases. The shock model also predicts 
that both CH + and H2O are rather well distributed across the var- 
ious velocity slices of the shock. This supports the assumption 
that the CH + and H2O blue line wings are similar, underlying the 
CH + line frequency determination. 

5.2. CH + absorption lines: turbulent dissipation models 

The large CH + column density observed in the foreground com- 
ponent corresponds to a los average of the CH + abundance, 
[CH + ]/[H] flve > 1.2 x 1(T 8 , for JV(H) = 1.4 X 10 22 cnT 2 
(see Sect. 2). It can be reproduced in the framework of TDR 
models. At any time, a very large number of tiny regions 
(I ~ 100 AU), altogether filling a small fraction of the en- 
tire los, are developing a transient warm chemistry triggered 
by dissipation bursts where local CH + abundances reach 10~ 6 
dGodard e t al. 2009). The filling factor of the ensemble of these 
tiny structures is set by the energy transfer rate in the turbu- 
lent cascade, e - pv 3 //, identified with the turbulent dissipa- 
tion rate. The resulting average abundance is found to scale 
as [CH + ]/[H] fll , c =6.4xl0- 9 (e/e )(« H /50cnr 3 r 2 - 6 for an ambi- 
ent radiation field x — 3 in ISRF units, where «h is the density 
of the gas in which the bursts occur and e is a non- local quantity 
of galactic average eo - 2x 10~ 25 ergcirT 3 s -1 (see lGodard et al.l 
2009). Fluctuations by two orders of magnitude are observed in 
the ISM about that galactic average, the highest values being 
reached in active star-forming regions, such as Cygnus X. We 
thus expect e to be up to 100 times higher than average in Cygnus 
X, so that the observed CH + abundance can be produced by in- 
termittent dissipation of turbulence occurring in gas of density 
up to ~ 250 cm -3 , The total number of tiny structures along that 
1.7 kpc los is found to scale as ~ 3.5 x 10 3 (e/eo)(«H/50cirT 3 )~ 2 . 



6. Conclusion 

We have detected the line profile of CFT(l-O) towards the DR 21 
massive star-forming region with the HIFI instrument, obtain- 
ing several major results. The line is a combination of broad 
emission and almost saturated broad absorption, which have 
comparable integrated areas. The rest-frame frequency of the 
CH + 7=l-0 line is inferred to be v = 835137 + 3 MHz. For 
the gas seen in emission, the CH + column density, of about 10 13 
cm -2 , compares well to predictions of C-shock models propagat- 
ing in dense, highly illuminated gas. Additional chemical mod- 
elling (and radiative transfer calculations) are needed to confirm 
the C-shock framework and exploit the other line signatures. The 
large CH + column density in the foreground gas, N{CH + ) ~ 10 14 
cm -2 , can be explained by the turbulent dissipation models in 
diffuse gas and confirms the large opacity of that line in the dif- 
fuse molecular gas, 
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